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We present 115 In NMR measurements in a novel thermodynamic phase of CeCoIns in high mag- 
netic field, where exotic superconductivity coexists with the incommensurate spin-density wave 
order. We show that the NMR spectra in this phase provide direct evidence for the emergence of 
the spatially distributed normal quasiparticle regions. The quantitative analysis for the field evo- 
lution of the paramagnetic magnetization and newly-emerged low-energy quasiparticle density of 
states is consistent with the nodal plane formation, which is characterized by an order parameter 
in the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state. The NMR spectra also suggest that the 
spatially uniform spin-density wave is induced in the FFLO phase. 



The interplay between magnetism and unconventional 
superconductivity with a nontrivial Cooper pairing has 
been a topic of recent intense study. The quasi-two 
dimensional heavy fermion superconductor CeCoIns [l[ 
continues to excite great interest, because it shows a 
number of fascinating superconducting properties 0-0|. 
Its superconductivity at high fields is destroyed by Pauli 
paramagnetic effect, as evidenced by the the first-order 
phase transition at the upper critical field H C 2 0-HI and 
anomalous flux line lattice form factor 0]. What is 
striking is that CeCoIns exhibits a new thermodynamic 
phase transition at (T*,H*) just below H c2 (Fig.QJa)) 
for both H \\ ab and H \\ c @, 0|. Closely related to 
the Pauli limited superconductivity, this high-field and 
low-temperature (HL) phase has been attributed to the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state 0, [Iq|, 
in which the pair-breaking arising from the Pauli effect is 
reduced by forming a new pairing state (k^,— k+ql) with 
nonzero q between the Zeeman splitted parts of the Fermi 
surface. One of the most fascinating aspects of the FFLO 
state is that Cooper pairs with finite center-of-mass mo- 
menta hq develop an oscillating superconducting order 
parameter in real space such as A(r) oc sm(q-r) and, as 
a result, nodal planes appear periodically perpendicular 
to the applied field [11] . 

The presence of the FFLO state in CeCoIns has been 
supported by several experiments, inclu ding impurity 
[l2|, pressure [T^j, and ultrasound studies fl^TflBI. How- 
ever, recent NMR fjB^ and neutron [13, Efl data 
in parallel field demonstrated a long-range static mag- 
netic order in the HL phase. The magnetic moment 
at Ce atoms is given by fi(r) = fiocos(Qs-r) with 
Qs = 2tt(£,£,^) (6 ~ 0.45) [H|, which is directed 
to the c axis, /jLq || c. Remarkably this incommensurate 
spin-density wave (IC-SDW) order vanishes when the su- 
perconductivity dies at H C 2 [HI, [l9j , indicating that the 
magnetism and the exotic superconductivity are closely 
intertwined. 

The observation of the magnetic order in the HL phase 
calls for a reexamination of a simple FFLO interpreta- 
tion. To account for the coexistence of d-wave supercon- 
ductivity and magnetic order confined exclusively in the 



superconducting state, several exotic superconducting or- 
ders, including pair-density wave state with a 7r-triplet 
component |2ll-l23j and spatially inhomogeneous SDW 
state induced around the FFLO nodal planes [24], have 
been proposed. However, the nature of the HL phase is 
still unclear and under hot debate. 

In this study, to improve our understanding of how 
the exotic superconductivity and magnetism can inter- 
act in CeCoIns, we measured NMR spectra on the three 
distinct In sites in parallel fields. The field evolution of 
the magnetization and the density of states (DOS) asso- 
ciated with the paramagnetic quasiparticle formation in 
the HL phase is extracted from the NMR spectra. Our 
results provide strong evidence for the formation of the 
FFLO state via the second order phase transition, which 
coexists with the static magnetic order. 

The 115 In NMR measurements were performed for 
H || [100] in the field-cooling condition by using a phase- 
coherent spectrometer on high-quality single crystals 
(T c = 2.3 K), whose physical properties are well char- 
acterized, as reported in the measurements of transport 
properties [25|, specific heat, magnetic susceptibility |j|, 
thermal conductivity [2] and ultrasound velocities [l4|. 
NMR spectra were obtained from a convolution of Fourier 
transform signals of the spin echo which was measured at 
40^50 kHz intervals. The Knight shift was obtained from 
central or satellite 115 In lines (J = 9/2) using a gyromag- 
netic ratio of 9.3295 MHz/T and by taking into account 
parameters of the nuclear quadrupole interaction. 

The tetragonal crystal structure of CeCoIn 5 consists 
of alternating layers of Celn3 and C0L22 (Fig.[]Jb)). In 
a magnetic field H || [100], there are three inequivalent 
In sites. The axially-symmetric In(l) is located in the 
Celn layer, whereas In(2a) and In(2b) sites are located 
between Co and Celn layers with the largest principal 
axis of electric field gradient parallel and perpendicu- 
lar to the applied field, respectively (Fig.QJc)). The 
NMR frequency in the superconducting state is spatially- 
distributed and is determined by the local magnetic 
field and hyperfine coupling to the conduction electron 
spins: H eS (r)=H+M (r)+A h fM s (r), where M (r) 
is the local magnetization due to the orbital (diamagnetic 



2 




FIG. 1. (Color online) (a) H-T phase diagram of CeCoIns 
for H || ab. (b) Crystal structure of CeCoIns. (c) A top 
view of the In (2) plane, and directions of the largest principal 
axis of electric field gradient (blue arrows) and of hyperfine 
fields (red arrows) on the In(2a) (cyan circles) and the In(2b) 
(blue circles) sites transferred from the magnetic Ce moments 
parallel to the c axis, (d) Field evolution of the NMR spectra 
at the In (2b) site as a function of K sp i n at T = 0.05 K in the 
normal (black), HL (red), and BCS (blue) states. i^ sp in = 
K — Kovh is estimated by using K G vh — 2.1%. Inset: Field 
dependence of the internal field H int = Sf /2^n obtained by 
the difference of the resonance frequency Sf of the two peaks 
(double-headed arrow in the main panel). 



screening current) effect, is the hyperfine coupling 
constant [26j] and Ms(r) is the local spin magnetization. 
The Knight shift, given as K = A^Ms/H, consists of the 
spin and orbital contributions, K = K sp { n + i\~ rt>- The 
amplitude of i\~ rb is estimated at low field and low tem- 
perature in the superconducting state where K sp [ n van- 
ishes. We note that in the vortex state of CeCoIns, line 
shapes at In(2a) and In(2b) sites with large amplitude 
of Ahf are predominantly determined by the hyperfine 
coupling (A M M s > M ). 

The antiferromagnetic staggered Ce moments due to 
IC-SDW with /xo || c induce the in-plane hyperfine field 
perpendicular and parallel to the applied field H at the 
In(2a) and In(2b) sites, respectively (Fig.QJc)), through 
dipolar type transferred hyperfine couplings [17j. Con- 
sequently, the In(2b) spectra should be shifted (into two 
peaks) by the antiferromagnetic staggered field, whereas 
In(2a) spectra should not. FigureQJd) depicts the field 
evolution of In(2b) spectra at 50 mK as a function of 
Knight shift. In the normal state for fi^H > 11.8 T, very 
sharp spectra with line width less than 50 kHz (black 
lines) are observed. For ji^H < 11.8 T the sharp In(2b) 
spectra broaden and split into two peaks with finite signal 
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FIG. 2. (Color online) Field evolution of the NMR spectra 
at the In(2a) site at T = 0.05 K in the normal (black), HL 
(red), and BCS (blue) states. The integrated intensity of each 
spectrum below H C 2 is normalized. The spin susceptibility 
(lower scale) is obtained by Xspm = i^spin/^hf with K OT h = 

I. 95%. The green shaded region indicates the quasiparticle 
spectrum emerged in the HL phase. Inset: Blow-up of spectra 
near the edge structure at fioH = 9.5 (black), 10.2 (blue), 

II. 1 (green), and 11.2 T (red). Dotted lines indicate the peak 
position in the normal state. 

weight between them (red lines), which is characteristic 
of IC-SDW long-range order along one spatial dimen- 
sion. In the BCS phase below H*(~ 10 T), NMR spectra 
(blue lines) becomes a single asymmetric line as expected 
in the usual vortex state. The inset of Fig.[TJd) shows 
the field dependence of the internal field iii nt determined 
by the difference of the resonance frequency of the two 
peaks in the HL phase. Above if*, H[ nt increases with H 
and jumps to zero at the first-order H C 2 transition after 
reaching maximum value of ~ 0.16 T, which corresponds 
to ~ 0.15/i#/Ce. The vanishing of SDW order above 
H C 2 and the magnitude of ifi nt are consistent with the 
previous NMR [l6|-[l8[ and neutron results [H,[20l. 

The field evolution of the In(2a) spectra is depicted 
in Fig.[2j The sharp resonance line in the normal state 
(black) becomes antisymmetric broad lines in the HL 
(red) and BCS (blue) phases. The most salient feature of 
the HL-phase spectra is the emergence of the edge struc- 
ture whose position coincides with the Knight shift in the 
normal state, as shown by the dotted line (also Fig. [21 in- 
set). This provides direct evidence for the emergence of 
the normal quasiparticle region in the HL phase |27l |28| . 
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FIG. 3. (Color online) Temperature evolution of the In(2a) 
spectrum in HL (red) and BCS (black) phases at /jloH = 
11.1 T. The peak intensity of each spectrum is normalized. 
Thin black lines (yellow shaded region) indicate the spec- 
trum at T — 0.21 K just above T*. The green shaded re- 
gion in the low-temperature data indicates the quasiparticle 
spectrum formed in the HL phase. 



We note that as the edge position is independent of mag- 
netic field, the effect of the hyperfine field due to the 
SDW ordering is negligibly small. This indicates that 
In(2a) spectra are predominantly affected by the local 
spin susceptibility arising from the normal quasiparticles. 

Here we comment on the reproducibility of the data. 
The NMR spectra are reproducible for different rf pow- 
ers and different cooling rates. We also measured the 
spectra on several different single crystals. Although the 
clear edge structure of In(2a) spectra is observed in all 
crystals, its sharpness slightly depends on the crystal. 
This appears to be related to the recent result that the 
HL phase is extremely sensitive to the nonmagnetic im- 
purity [12| . The NMR spectra with less pronounced edge 
structure can also be found above H* in Ref . fl8l. 

We stress that the clear-cut sharp edge structure is 
consistent with the recent calculation of the NMR spectra 
in the FFLO state based on the microscopic Bogoliubov- 
deGennes equation [28|. For the quantitative analysis, 
we measured the temperature dependence of the In(2a) 
spectrum (Fig. [3]). As shown by thin black line, which 
is the spectrum at T = 0.21 K just above T*, the shape 
of the main peak remains almost unchanged with tem- 
perature. This allows us to separate the spectrum in the 
edge region (the green shaded region in Fig. [3]) by sub- 
tracting the BCS spectrum at 0.21 K (the yellow shaded 
region) from the whole spectra in the HL phase. The 
green shaded regions in Fig. [2] also indicate the spectra 
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FIG. 4. (Color online) (a) Field dependence of the param- 
agnetic magnetization M p of the In (2a) (red circles) and the 
In(l) (blue open circles) sites. Inset: Schematic quasiparticle 
structure in the FFLO state. The nodal planes (green) with 
period of 2ix /q appear perpendicular to the Abrikosov vortex 
lattice (blue), (b) Field dependence of the DOS of the param- 
agnetic quasiparticles, n qp and In, extracted from the In(2a) 
spectra. The solid line is a fit to \/H — independence. 



near the edge regions obtained by the similar method. 

The total paramagnetic magnetization M p 
is evaluated from the whole NMR spectrum 
P(-Kspin) by integrating the spin part of the 
Knight shift as M P (H) = -*L g (H), where 
g(H) = J K spin P(K spin )dK spin / fp(K spin )dK spin . 
FigurelUa) depicts the field dependence of M p obtained 
from the In(2a) spectra as well as M p (H) obtained from 
the same analysis for the In(l) spectra. Both M p (H) 
from the In(l) and In(2a) spectra increase linearly with 
H in the HL phase. The difference of the paramag- 
netic magnetization between the HL and BCS phases, 
SM P (H) = M P (H) - M£(H), increases continuously 
from zero as SM p (H) oc (H — H*) in the HL phase, 
where M^{H) is obtained by a linear extrapolation 
from the BCS phase (dotted line in Fig. Ufa)). We 
emphasize that M P (H) cannot be obtained from the 
bulk magnetization measurements [4j which contains the 
contribution of the antiferromagnetic Ce moments. 

The DOS of the normal quasiparticles emerged in the 
HL phase is proportional to the area of green shaded re- 
gions in Figs. [2] and [31 as the intensity there should be 
proportional to the number of nuclei which detect the 
quasiparticle susceptibility. FigureH^b) depicts the field 
dependence of the normalized DOS, n qpi which is ob- 
tained by the area of the green region divided by the 
area of the full spectrum at each field in Fig. [2l To es- 
timate the DOS in an alternative way, the intensity In 
of the spectrum at the normal state Knight shift (dotted 
line in Fig. [2]) is also plotted in Fig.|4jb). Both n qp and 
In increase in proportion to \JH — H* in the HL phase. 

The i^-linear dependence of the paramagnetic magne- 
tization and V^-dependence of the low-energy DOS both 
provide key information about the order parameter that 
characterizes the emergence of the normal quasiparticles 
in the HL phase. We stress that both field dependencies 
are exactly what are expected in the FFLO state where 
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the order parameter is described by the amplitude of the 
modulation wave vector q =\q\. The change in the para- 
magnetic magnetization, which is the first derivative of 
the free energy with respect to the applied magnetic field, 
is proportional to the square of the order parameter near 
the second order phase transition. Then it increases in 
proportion to the magnetic field as SM P oc q 2 oc (H — H*) 
near the FFLO transition [H, Hq|. Moreover, the DOS 
of the normal quasiparticles emerged in the FFLO phase 
is proportional to the number of nodal planes, which is 
proportional to q (inset of Fig.H^a)); n qp is expected to 
increase with H as n qp oc q oc \JH H*. Thus the quan- 
titative analysis of the In(2a) and In(l) spectra provide 
strong support for the formation of the FFLO state. 

It has been theoretically proposed that the IC-SDW 
moment is induced around the FFLO nodal planes 
In the presence of such a spatially non-uniform IC-SDW 
state, the NMR spectra at the In(2b) site in the HL phase 
consist of both contributions of the regions far outside 
(BCS spectrum with one peak) and around the nodal 
planes (IC-SDW spectrum with two peaks). In contrast, 
the present In(2b) spectra in the HL-phase shown in 
Fig.QJd) consist entirely of the split one with no dis- 
cernible component of the BCS spectrum. This suggests 
rather uniform IC-SDW, in which the magnetic order is 
present even in the region far away from the nodal planes. 
This is consistent with the absence of satellite peaks in 
the neutron scattering experiments [20j . 

We point out that the present results are incompatible 



with the pair-density wave scenario [2l|, [22| . The emer- 
gence of the normal quasiparticles in the HL-phase is not 
expected in this scenario. Moreover, the order parameter 
g, which depends on H in accord with the FFLO second 
order transition, is at odds with this scenario that as- 
sumes the Cooper pair with the SDW modulation wave 
vector Qs, which is field independent [19j . 

It has been reported that the magnetic moment in- 
duced by the SDW disappears when the magnetic field 
is tilted away from the ab plane by 17° [3l|. Moreover, 
a possible FFLO phase with no magnetic order appear 
in H || c d, These results appear to indicate that 
coupling between FFLO and IC-SDW becomes weaker 
with tilting H from the ab plane. Recently, it has been 
suggested that in CeCoIns with ^2 -wave symmetry 
the strong Pauli paramagnetism plays an important role 
for the SDW formation as well as the FFLO state par- 
ticularly in parallel field [33|, [34| . Further investigation 
is required to understand the perplexing relationship be- 
tween the coexisting FFLO and IC-SDW states. 

In summary, 115 In NMR measurements demonstrate 
the emergence of a spatially distributed normal quasipar- 
ticle region in the HL phase of CeCoIns in parallel field. 
The field evolution of the paramagnetic magnetization 
and low-energy quasiparticle DOS can be described well 
by the order parameter associated with the nodal plane 
formation via the FFLO second order phase transition. 
The NMR spectra also reveal that the spatially uniform 
SDW coexists with the FFLO nodal planes. 

We thank N. Curro, R. Ikeda, K. Ishida, K. Machida, 
H. Shimahara, and Y. Yanase for valuable discussion. 
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